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Abstract 

Recent results from lattice QCD simulations provide a realistic pic- 
ture, based upon first principles, of T physics. We combine these results 
with the experimentally measured mass of the T meson to obtain an ac- 
curate and reliable value for the fo-quark's pole mass. We use two different 
methods, each of which yields a mass consistent with M& = 5.0(2) GeV. 
This corresponds to a bare mass of M° = 4.0(1) GeV in our lattice theory 
and an MS mass of Af 6 MS (M 6 ) = 4.0(1) GeV. We discuss the implications 
of this result for the c-quark mass. 

PACS numbers: 12.38. Gc, 14.40. Gx, 14.56.Fy, 12.38.Bx 

The mass Mb of the b quark is a fundamental parameter of the Standard 
Model. Various phenomenological studies have suggested masses in the range 
of 4.5-5 GeV [Q. In this paper we use new results from lattice QCD to obtain a 
value for Mi, from the measured mass of the T meson |2| . We employ two meth- 
ods for computing the quark mass; both are consistent with M b = 5.0(2) GeV, 
which is presently among the most accurate and reliable of determinations. 

The notion of quark mass is complicated by quark confinement ||. In a 
particular process, a 6-quark's effective, or running, mass Mb(q) depends upon 
the typical momentum q transferred to it. In perturbation theory, at least, this 
mass stops running when q falls below Mb(q). Below this point, Mf,(q) equals 
the quark's pole mass, which is defined perturbatively, in terms of the bare mass, 
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by the pole in the quark propagator. It is this mass we quote for Mb above. The 
use of perturbation theory is certainly justified when q = Mb. At very low q's the 
effective mass may be modified away from our value by nonperturbative effects, 
but our value should hold for Aqcd Cg<Mj. It is therefore the appropriate 
mass to use in quark potential models and other phenomenological applications 
that involve this momentum range. We may also use perturbation theory to 
relate our bare quark mass to the MS mass; we find M b MS (Alb) — 4.0(1) GeV. 

Both of our mass determinations rely on numerical simulations of the T using 
lattice QCD. In lattice QCD spacetime is approximated by a discrete grid of 
finite physical volume. An action for QCD is defined in terms of quark and gluon 
fields on the nodes of the grid and the links joining them. Monte Carlo methods 
are used to evaluate the path integrals that define various correlation, or Green's, 
functions. The energies and wavefunctions of the T and its excitations are 
extracted from these correlation functions. 

Our gluon configurations were generated using the standard Wilson lagrangian 
for gluonsHl. For most of our simulations we used a 16 3 x 24 grid with j3 = 
6/g 2 = 6 ||, where g is the bare coupling constant. This corresponds to an 
inverse lattice spacing of a -1 = 2.4(1) GeV and a volume of approximately 
1.3 3 x 2.0 fm 4 . The lattice volume is more than adequate since T's have a ra- 
dius of only i?x ~ 0.2 fm. The finite grid spacing introduces systematic errors 
of order a 2 / R\] we have estimated these to be less than 5% (for binding ener- 
gies) both from quark potential models and from other simulations with smaller 
lattice spacings ||. 

The gluonic configurations described above do not include effects due to 
light-quark vacuum polarization. There are strong theoretical arguments sug- 
gesting that T physics is largely unaffected by the light quarks, at least for states 
that are well below the BB threshold. To verify this we repeated a part of our 
analysis for gluon configurations that contain contributions from n/ = 2 flavors 
of light quark |jj . The mass of the light quarks in these configurations (under 
100 MeV) is larger than the mass of a u or d quark but negligible relative to 
the typical momentum transfers in an T (about 1 GeV). Thus the quarks are 
effectively massless. The lattice spacing and volume are the same for these con- 
figurations as for the n/ = configurations described above. The bulk of the 
results presented in this paper are based upon the simulations with nf — 0; how- 
ever, we found these essentially unchanged by the inclusion of light-quarks 

We used the NRQCD lagrangian for the b quarks This is a nonrela- 
tivistic lagrangian that is tuned to efficiently reproduce results from continuum 
relativistic QCD order-by-order in the lattice spacing a and quark velocity v. 
We included the leading relativistic and finite-lattice-spacing corrections; further 
corrections are almost certainly negligible for our purposes. The only parameter 
in this lagrangian, other than the QCD charge g, is the bare quark mass M t °. 
This is tuned so that the simulation gives the correct dispersion relation for 
the T. The energy of a low-momentum T can be computed in the simulation, 
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and has the form 



Ey(p) « -Bnr(T) + p 2 / 2M kin (T) , (1) 

where £/nr(T) is related to the nonrelativistic binding energy of the meson, and 
•Mkin(T) is the kinetic mass of the meson. We tuned the bare quark mass M" 
until Mkin(T) was equal to the measured mass of the T.[] We found that the 
correct bare mass is given by oM° = 1.7(1). 

Our simulations very successfully reproduce the general features of T physics. 
This is illustrated by Fig. [I], which compares simulation results for the spectrum 
of the lowest-lying states with experiment. In Fig. ^, simulation results for the 
fine structure of the lowest lying P state are compared with experiment. The 
excellent agreement supports the reliability of these simulations. We emphasize 
that these are calculations from first principles; our approximations can be 
systematically improved. The only inputs are the lagrangians describing gluons 
and quarks, and the only parameters are the bare coupling constant and quark 
mass. In particular, these simulations are not based on a phenomenological 
quark potential model. 

Our first procedure for computing the 6-quark mass uses simulation results 
for the nonrelativistic energy Snr(T) of the T (Eq. ([!])). The quark mass is 

M b = (1/2) (M T - (E NR - 2Eo)) , (2) 

where Mx =9.46 GeV is the experimentally measured mass of the T, and Eq 
is the nonrelativistic energy of a & quark with p = in NRQCD. The quantity 
-E'Nr(T) — 2Eq is the effective binding energy of the meson. The quark energy Eq 
is ultraviolet divergent as the lattice spacing a vanishes, and so at small a, it 
depends primarily on physics at momenta of order l/a. Thus it can be computed 
using weak-coupling perturbation theory pi, with 

£ = (b Q /a) a v {q ) (1 + 0{a v )) ■ (3) 

Here ay(qo) is the strong coupling constant defined in |p"i"f j 2 ] 

In Table @ we present our simulation results for E'nr, and the corresponding 
.Eo's from perturbation theory. We list results for several values of the bare 
quark mass; for the 6-quark, aM® = 1.7(1). From these results we conclude 
that the pole mass of the b quark is Mb = 5.0(2) GeV. The major sources of 
uncertainty in this determination are: 

'In a purely nonrelativistic theory, Galilean invariance implies that the kinetic mass of 
a bound state equals the sum of the masses of its constituents. In a relativistic theory the 
binding energy enters as well, resulting in a 10% shift for T's. Thus it is important to include 
relativistic corrections in the quark action when tuning the quark mass using Afkin. 

jiAU perturbative results in this paper are obtained using the techniques and conventions 
of (ill ; in particular, the definition of the running coupling constant and the procedure for 
setting the coupling scale are from this reference. The value of ay is determined from the 
measured expectation value of the plaquette operator, as discussed in |ll| . This coupling 
constant is related to the MS coupling by: otr=-(Q) = »^(exp(5/6) Q) (1 + 2aty/n). 
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Figure 1: NRQCD simulation results for the spectrum of the T system including 
radial excitations. Experimental values (dashed lines) are indicated for the 
triplet 5-states, and for the spin-average of the triplet P-states. The energy 
zero from simulation results is adjusted to give the correct mass to the T(l 3 S'i). 
These data are for n/ = 0. 
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Figure 2: Simulation results for the spin splittings between the lowest lying 
P-wave states in the T family. The dashed lines are the experimental values for 
the triplet states. Energies are measured relative to the center of mass of the 
triplet states. These data are for n/ = 0. 
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n f 


aM° 


aE NR (T) 
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aE a 
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6.0 





1.71 


0.453(1) 


0.75 


0.48 


0.32 


4.95 






1.80 


0.451(1) 


0.77 


0.51 


0.31 


4.92 






2.00 


0.444(1) 


0.81 


0.55 


0.30 


4.91 






3.00 


0.412(1) 


0.91 


0.68 


0.28 


4.92 


5.6 
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1.80 


0.493(1) 


0.77 


0.51 


0.38 


5.06 


6.0 





1.80 


0.448(1) 


0.76 


0.47 


0.33 


4.98 



Table 1: Simulation and perturbative results used in the first method for de- 
termining Mb. Values for Mb are in GeV, and are determined using Eq. (6). 
Results are presented for n/ = and 2 light-quark flavors, and for a range of 
bare quark masses; the 6-quark has aM® = 1.7(1). Results are also given for 
just the leading term in the NRQCD lagrangian (SH = 0). 



(1) Tuning errors in the bare quark mass. The binding energy, -Enr — 2i?o, 
like the radial and orbital level splittings, should be almost independent of the 
quark's mass for masses of order Mb. Thus Mb, as determined from Eq. (||), 
should also be independent of the quark's mass. This is confirmed by our results 
in Table 13. Errors in Mb due to uncertainty in the bare quark mass are less 
than 1%. 

(2) Two-loop corrections to Eq. Quantitative studies jll] of the reliability 
of perturbation theory at scales of order qq suggest that two-loop corrections 
to Eq could range from 1-30% of Eq. A 20% contribution would shift Mb by 
about 3%. Two-loop corrections could not be much larger than this without 
ruining the agreement we find between results obtained with widely different 
quark masses. Nonperturbative effects are most likely smaller. 

(3) Simulation errors in £'nr(T). The dominant sources of error in -Enr. are 
the 0(a 2 ) errors in the gluon action, and our use of zero and two flavors of light 
quarks rather than three. Based upon the results in the Table [j], together with 
results from simulations with smaller lattice spacings |(| , we expect these effects 
to shift the quark mass by less than 100 MeV. Statistical errors are negligible 
here. 

(4) Simulation errors in a~ 1 . The lattice spacing is determined by matching 
simulation results for the T spectrum to experiment. It is easily measured to 
within ±5%. Since -EfjR an d 2i?o are each only about 10% of Mt, this amounts 
to an uncertainty of less than 0.5% in Mb. 

The relativistic and finite-lattice-spacing corrections in the quark lagrangian 
should have only a small effect on the total mass of the T. This is confirmed by 
our simulation, as illustrated in Table 0. The last entry in the table was obtained 
from a simulation without these corrections. It agrees well with results from the 
full simulation, confirming that further corrections are probably unimportant. 
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Our second procedure for determining Mb is to tune the bare quark mass Af" 
until the kinetic mass of the T, as computed in the simulation (Eq. (|l|)), agrees 
with the T's true mass. Then the pole mass is Z m M®, where the renormaliza- 
tion constant Z m is computed using perturbation theory [ |To[ , with 

Z m = l + b m a v {q m ) + 0(a 2 v ). (4) 

To reduce the sensitivity of the result to the values of the lattice spacing and 
bare quark mass, we rewrite the expression for the pole mass as 

n M 

Here My = 9.46 GeV is the T's experimental mass, while Mk; n (T) is its 
mass as determined from a simulation using the bare mass M". Our results 
are summarized in Table ||. For n/ = 0, the kinetic mass is closest to the 
T mass when aM® = 1.71, and therefore the renormalized quark mass is 
Mb = 4.9(2) GeV. The corresponding bare quark mass is M° = 4.1(1) GeV 
(from Eq. (^|) with Z m — > 1). The nj = 2 data must be extrapolated slightly in 
aM® to obtain the correct T mass. Using the rif — data as a guide for this 
extrapolation, we obtain M b = 5.0(2) GeV and M° = 4.0(1) GeV for n f = 2. 
The main sources of uncertainty in these determinations are: 

(1) Two-loop corrections to Z m . The one- loop corrections to Z m shift Mb by 
15-20%. If two-loop corrections are 15-20% of the one-loop corrections, they 
could shift Mb by 2-4%. The bare mass is unaffected by these corrections. 

(2) Simulation errors in aMki n (T). These are mainly due to the omission of 
0(v 4 ) relativistic corrections in the quark action. Such effects should shift Mb 
by less than 1%. (The 0(v 2 ) corrections, which we include, shift Mb by only 
7%.) Statistical errors in aM\c iri are of order 1%. 

(3) Tuning errors in the bare quark mass. Since a -1 is known only to ±4%, 
the T's kinetic mass can only be tuned with an accuracy of ±4%. However our 
determination of Mb (Eq. (^|)) minimizes sensitivity to such tuning, as is evident 
from Table E| Thus tuning errors have a negligible effect on Mb- 

In this paper, we have presented two different determinations of the 6-quark 
mass, each consistent with a pole mass of Mb — 5.0(2) GeV. The systematic 
errors are quite different in each case, and have been probed by computing 
with a variety of quark masses and quark lagrangians, with and without light- 
quark vacuum polarization. The complete agreement between the two methods 
is a strong indication of the validity of each, and more generally of the lattice 
QCD techniques on which they rely. In both cases the dominant uncertainties 
are in the perturbative calculations. A two- loop evaluation of either the quark 
energy Eq or the mass renormalization constant Z m would probably reduce 
the errors in Mb by a factor of two or more. We are currently examining the 
feasibility of such calculations. 



7 






n f 


aM° 


aM kin (T) 




aq m 


Z m ~ 1 


Af kin (T) 


M b 


M° b 


6.0 





1.71 


3.94(3) 


0.48 


0.49 
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4.49(4) 


0.42 
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0.16 


10.7(4) 


4.88 


4.21 






3.00 


6.57(7) 


0.24 


0.52 


0.09 


15.8(6) 


4.73 


4.32 


5.6 


2 


1.80 


4.17(2) 


0.46 


0.51 


0.23 


10.0(4) 


5.02 


4.08 



Table 2: Simulation and perturbative results used in the second method for 
determining M b . Values for Mk; n (T), M b and M° (last column) are in GeV. 
Values for M b are determined using Eq. (9). Values for M 6 ° are obtained using 
Eq. (9) with Z m — > 1. Results are given for rif = and 2 light-quark flavors, 
and for a range of bare quark masses; the 6-quark has aM? = 1.7(1). 

Combining the standard result for the pole mass in terms of the MS mass Jl2[ , 

M b = M b ME {M b ) {l + 0.424 ay (0.22M b ) + 0.164 a 2 v + •■■}, (6) 

with our result in terms of the bare quark mass, M b = Z m M^, we obtain a 
relation between the bare quark mass on the lattice and the MS mass: 

Mp(M 6 ) = M° {1 + 0.06a v (q m ) + ■••}, (7) 

where w 0.8/a for our Af°'s. Using this relation, our data implies that 

M b MS (Af b ) = 4.0(1) GeV. 0(c? v ) corrections should be smaller than in Eq. (§) 
for the pole mass. This is because the scale for ay is significantly larger here 
and therefore ay is smaller. The larger scale is sensible since the MS mass is a 
bare mass and so should be less infrared sensitive than the on-shell pole mass. 

Our 6-quark mass can be used to estimate the c-quark's pole mass. Because 
of heavy-quark symmetry, the difference between the spin-averaged mass of 
the B and B* mesons, M B i B * , and M b is the same as the difference between 
the spin-averaged D/D* mass and M c up to corrections of order Aq CD /M c . 
Thus we expect 

M c w M D/D , + (M b - M B/B .) - 1.6(2) GeV. (8) 

A different estimate is possible if the binding energies of the T and ip are roughly 
equal. This is probably the case since radial and orbital excitation energies are 
roughly equal in the two systems. Also our simulation shows that the binding 
energy is almost mass independent for masses larger than M b . Thus we expect 

M c w (1/2) (Mip + (2M 6 - M T )) = 1.8(2) GeV. (9) 

We are currently extending our lattice analysis to allow a direct determination 
of M c . 
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